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We present a new method of gauge fixing to standard lattice Landau gauge, Max Re Tr ^2 x U^, x , in which the 
link configuration is recursively smeared; these smeared links are then gauge fixed by standard extremization. The 
resulting gauge transformation is simultaneously applied to the original links. Following this preconditioning, the 
links are gauge fixed again as usual. This method is free of Gribov copies, and we find that for physical parameters 
(j3 > 2 in SU(2)), it generally results in the gauge fixed configuration with the globally maximal trace. This method 
is a general technique for finding a unique minimum to global optimization problems. 



There are two outstanding difficulties with 
gauge fixing on the lattice: 

1. The Gribov problem: covariant gauge fixing 
inevitably leads to multiple solutions of the 
gauge fixing condition. 

2. The Smoothness condition: if one alters the 
gauge condition so that it is free of Gribov 
copies (like axial gauge) it is generally not 
smooth, and difficult to compare to pertur- 
bation theory. 

Below we present a method addressing these 
two issues which is both simple and fast. It was 
inspired by the use of improved operators and 
the smoothing properties of smearing links. It 
is free of Gribov copies, and usually produces 
the smoothest (Max I{U^G) = Re Tr J2 x fl U°) 
configuration. 

It should be contrasted to other covariant so- 
lutions to the gauge fixing problem [jjj, which 
are also Gribov copy free, but require non-trivial 
computational resources and do not rotate fields 
to Landau gauge. 

1. Method 

We begin by iteratively smearing a copy of the 
links, ie. recursively smearing again the smeared 
links, with an APE type smearing process. In 
this study we use the method of j^j in which links 
are averaged with connecting staples; staples have 



weight w, and the resulting link is reunitarized 
into SU(N). 

We have found numerically that as we iterate 
this process there is a critical weight w c below 
which the configuration is cooled, ie. the average 
plaquette, moves closer and closer to the iden- 
tity; with a weight above w c smearing moves the 
links to a "rough" configuration with Tr(n) ~ 0. 
The value of this critical weight seems to depend 
weakly on /3 and volume: w c ~ 0.5 — 0.6. 

This method of cooling however is equivariant: 
If starting configurations U A and U B are related 
by a gauge transformation G AB , then the corre- 
sponding smeared configurations and V? are 
related by the same gauge transformation G AB . 
This is useful since the smeared configuration is 
approaching the trivial orbit, where we know the 
unique Landau gauge configuration. It is the con- 
figuration in which each link is constant and diag- 
onal: U^x) = (Df,) 1 /^ such that TrD^ = TrP M 
where is the smeared Polyakov loop in the 
^-direction (made of n M links). We neglect here 
issues of degeneracy of eigenvalues, etc. 

Thus if we smear extremely close to the triv- 
ial orbit, we can find a Gvd which rotates the 
smeared links to Landau gauge uniquely, then 
apply this Gvd to the original links U^, produc- 
ing a unique configuration D'^; due to the non- 
linearity of non-abelian gauge transformations, it 
is neither diagonal nor satisfies the Landau gauge 
condition in SU(N) theory. Still it represents a 
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unique starting point on the physical orbit, from 
which we can gauge fix to Landau gauge and be 
assured of a unique solution. This is displayed in 
figure 1. 
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Figure 1. Gauge fix to £> M and use the same 
G to fix (7 M to D'p. 

2. A simpler way 

In our experiments we have found that it is 
not necessary to smear the links all the way to 
the trivial orbit; after only a moderate amount 
of smearing, gauge fixing of the smeared links is 
unique. Our view of this is captured below in an 
artist's rendering of the change in I(V,G) as we 
smear. 

If the original configuration is smooth enough, 
the global maximum on the smeared orbit may 
lie in the basin of attraction of the (gauge trans- 
formed) global maximum on the physical orbit. 
When this happens, the preconditioning step of 
gauge fixing the smeared links to Landau gauge 
leads to the final copy being the one of maxi- 
mal I(Ufj,,G), ie. the smoothest copy. If the 
original configuration is quite rough (unphysically 
so) , there often is a mismatch between the global 
maxima of the smeared and original / functions. 
The resulting configuration then does not have 
the globally maximal trace, although the trace is 
usually still rather high among the copies. In ei- 
ther event though, gauge fixing (to Landau 
gauge) is unique. 
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Figure 2. The gauge fixing functional — I(V, G) 
changes as we smear V^. Initially the func- 
tional has many relative minima leading to Gri- 
bov copies. As the links are recursively smeared 
— I(V, G) becomes smoother, developing a unique 
minimum. 

3. Parameters 

In the original form of this method, we should 
smear until we are within machine zero of the triv- 
ial orbit (F^ = 0), however we have found that 
only partial smearing is necessary. The optimal 
parameters specifying this have emerged rather 
hcuristically, and we will simply summarize them 
here and refer the reader to our more complete 
paper on this method of gauge fixing || . Also we 
have studied small lattices so far, up to 12 4 , and 
we gauge fix to variations smaller that 10~ 13 in 

The smearing coefficient w, which is the weight 
given to the staples should be as large as possible, 
but below the critical value so that the smeared 
configuration is smoothed. We have found that 
w ~ 0.3 works best, looking at how topological 
objects (which slow the cooling) are annihilated. 

The most important parameter is the number 
of smearing sweeps done. We have found that this 
depends on the number of topological objects in 
the configuration. Thus smearing is pursued until 
the action falls below the one-instanton bound. 
On our 8 4 lattices, we smeared until the average 
plaquette was within 10~ 5 of 1, which then leads 
to unique gauge fixing. In figure 3 we show the 
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cooling history of 6 lattices at a range of j3 values. 
The plateaux show the stability of instanton pairs 
under smearing (like under cooling). 
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Figure 3. 2.0 — Tr(d) versus smearing (w = 0.3) 
on lattices of size 8 4 . 



4. Summary 

In figure 4 we display typical results on 100 
random gauge transformations at 3 f3 values: 2.0 
(top), 1.75 (middle), and 1.5 (bottom). The 
graphs contain histograms of the final value of 
/ on the gauge fixed copies obtained using local 
cxtrcmization, and the arrow shows the unique 
copy obtained with our method. To summarize 
our method: 

• Recursively smear a configuration into 

• Gauge fix by extremizing, I(V, G) = 
£ Xj/J Re Tr GHx)V^(x)G(x + ft), and ap- 
ply the same gauge transformation to U^. 

• Finally gauge fix as usual, by extremiz- 
mg I (U,G). 

This construction trivially generalizes to 
Coulomb or maximal Abclian gauges, and can 
be more generally applied to some optimization 
problems. 



0.7848 0.785 0.7852 

beta - 1 .75 



3.7854 0.7S56 





Figure 4. Gribov copies obtained by standard ex- 
tremization and the unique copy j obtained from 
smearing. 
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